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Abstract 

Using data collected by the CLEO detector in the T(45) region, we report 
new measurements of the exclusive decays of B mesons into final states of 
the type Afpn(n), where n = 0, 1,2,3. We find signals in modes with one, 
two and three pions and an upper limit for the two body decay A+p. We 
also make the first measurements of exclusive decays of B mesons to S c pn(7r), 
where n = 0, 1, 2. We find signals in modes with one and two pions and an 
upper limit for the two body decay S c p. Measurements of these modes shed 
light on the mechanisms involved in B decays to baryons. 
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A distinctive feature of the B system is that the large mass of the 6-quark allows weak 
decays of the B mesons to proceed via the creation of a baryon anti-baryon pair. The mech- 
anisms for baryon production have been the subject of several studies in the last decade. 
The dominant decay mechanism is expected to be b — > cud transitions via internal or ex- 
ternal W^-decay. This can lead to final states including a charm meson as well as a baryon- 
antibaryon pair, as recently observed by CLEO |l] and Belle 0. However, the simplest 
decay diagrams lead to states of the form B —> A C NX, where N represents an anti-nucleon. 
Charge-conjugate processes are implied throughout this paper. Inclusive studies of A+ pro- 
duction from B decays indicate a branching fraction of around 5%, and the soft A+ mo- 
mentum spectrum indicates that multi-body decays dominate In 1996, CLEO made 
the first exclusive measurements || of decays of this type, and found B(B~ — > A+p7r~) = 
(0.62±g;|? ± 0.11 ± 0.10) x 10~ 3 and B(B° -> A+pvr+Tr-) = (1.33^42 ± 0.31 ± 0.21) x 10~ 3 . 
The analysis presented here uses a larger CLEO data sample and improved analysis tech- 
niques to make further measurements of this type, confirming the previous observations and 
measuring new modes B~ — ► A^pn + n^n~ and B~ — > A+p7r~7r°. Furthermore, by investi- 
gation of the resonant substructure of these decays, the first exclusive decays of the form 
B — > H C NX have been measured. Comparisons of the branching fractions of these modes 
gives information on the underlying mechanisms involved. 

The data were collected with two detector configurations, CLEO II g and CLEO II. V 
at the Cornell Electron Storage Ring. The data comprise 9.17 fb^ 1 taken at the T(4S) which 
corresponds to 9.74 x 10 Q BB pairs, together with 4.6 fb~ taken in the e + e~ continuum below 
the T(AS) that are used to evaluate possible backgrounds. We assume that the produced 
B + B- rate is the same as B°B~® at the T(4S r ). 

The signal B meson candidates are fully reconstructed by combining detected photons, 
protons, and charged kaons and pions. The tracking system consisted of several concentric 
detectors operating inside a 1.5 T solenoid. For CLEO II, the tracking system consisted of 
a 6-layer straw tube chamber, a 10-layer precision drift chamber, and a 51-layer main drift 
chamber. The main drift chamber also provided a measurement of the specific ionization, 
dE/dx, used for particle identification. For CLEO II. V, the straw tube chamber was replaced 
by a 3-layer, double-sided silicon vertex detector, and the gas in the main drift chamber was 
changed from argon-ethane to helium-propane mixture. In both configurations, photons were 
detected by an electromagnetic calorimeter consisting of 7800 Thallium-doped Csl crystals. 

A+ candidates are reconstructed in the modes pK~n + ,pK s and Air + , where K s — > 7r + 7r~ 
and A — ► pir~. The Kg and A candidates are reconstructed from oppositely charged tracks 
which form a vertex well detached from the main event vertex in the plane transverse to 
the beam direction. The invariant mass of the Kg(A) candidate is required to be within 8 
(3.5) MeV/c 2 of the known mass. Neutral pion candidates are formed from pairs of showers 
detected in the calorimeter which yield a 77 invariant mass within 2 standard deviations of 
the known 7r° mass. The A, Kg and 7r° candidates were then all kinematically constrained 
to their known masses. 

Particle identification of p, K~, and ir + candidates was performed using specific ionization 
measurements in the drift chamber, and when present, time-of-flight measurements. For each 
mass hypothesis, a combined % 2 probability Pj was formed (i = ir, K,p). Using these Pj's, a 
normalized probability ratio Li was evaluated, where Li — Pi/ (P n + Pk + P p ). Real protons 
have L p of close to 1.0, whereas tracks due to other particles are concentrated near L p = 0.0. 
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For a track to be used as a proton daughter of a A+ —» pK~rr + or pK$ decay, we require it to 
have L p > 0.9, which eliminates much of the background but with considerable diminution of 
efficiency. For kaons we applied a slightly looser and more efficient requirement of Lk > 0.7. 
We have chosen these selection criteria using a Monte Carlo simulation program to maximize 
the significance of the A+ signals. The candidate anti-proton which is the direct daughter of 
the B meson has a looser requirement of L p > 0.1 as it has a higher momentum distribution 
and lower backgrounds than the decay daughter of the A+. The proton from the A — > piT~ , 
and all the pion candidates, are required to have particle identification parameters consistent 
with their hypothesis. Tracks with no particle identification information are assumed to be 
due to pions. 

To suppress continuum background, the normalized Fox- Wolfram second moment |§ is 
required to be less than 0.35. The number of A+ candidates from the T(45) data, after the 
contribution from continuum events is accounted for by a subtraction of scaled continuum 
data, is 7100 ± 350 (12100 ± 450) from the CLEO II (CLEO II. V) detector configurations. 

To reconstruct exclusive B decays we select A+ candidates whose mass is within 2a of the 
nominal mass. The mass resolution, a, was calculated for each of the three decay modes and 
two detector configurations separately by use of a GEANT-based Monte Carlo simulation 
program |§. We constrain the mass of these A+ candidates to the A+ peak value using a 
kinematic fitting program, and combine them with an anti-proton candidate and a number 
of pion candidates. We define the beam-constrained mass as Mb = \J E~1 carn — (Z^p.;) 2 , where 
Pi is the 3-momentum vector for the i th daughter of the B candidate and -Ebcam is the beam 
energy. The resolution in Mb is dominated by the spread in the CESR beam energy and is 
much better than the resolution in the invariant mass of the combination. 

For each combination, we calculated the energy difference AE = E meas — E heam , where 
-E'meas is the measured energy of the combination. A correctly reconstructed B meson has a 
AE distribution with a maximum at GeV. The AE resolution, <jae, was calculated for 
each mode and detector configuration separately using the Monte Carlo simulation program, 
and combinations were required to have \AE\ < 2<Jae- A further reduction in background is 
achieved by cutting on Qg, the polar angle of the B in the laboratory frame with respect to 
the e + e~ axis. The distribution of cosO^ is proportional to sm 2 Q B for e + e _ — > T(4S f ) — > BB, 
whereas background events are distributed nearly isotropically. We require |cos6b| < 0.9. If 
there are multiple candidates in an event with M B > 5.2 GeV/c 2 for a given decay channel, 
the entry with the smallest absolute value of AE is selected. 

The Mb distributions, after all selection criteria have been applied, are displayed in Fig. [I] 
for all modes investigated. Strong signals are found in the modes A+p7r~ and A+p7r~7r + , 
confirming the previous observation of these modes. Signals are also found in the new modes 
A+p7r~7r + 7r~, which has a 13.0a significance^ and A+p7r _ 7r°, which has a 8.2cr significance. 
There is no statistically significant signal in the two-body decay B° — > A+p. Each Mb 
distribution is fit to a fixed width Gaussian signal function, and a background function of 



1 We define our significance as the probability, expressed in normal distribution sigma, of our 
expected background to fluctuate to our signal's central value. Poisson (Gaussian) statistics are 
used for expected backgrounds with less (greater) than 30 events. 
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an exponential with phase-space threshold suppression. The signal yields from these fits are 
shown in Table |, where the uncertainties are statistical only. We have verified that similar 
distributions made with A+ sidebands, AE sidebands, or continuum data show no peaking 
in the B mass region. 

Knowledge of the substructure of the multi-particle final states is very important. From 
a purely practical point of view, the substructure changes the efficiency for detecting a final 
state. This is particularly true when the intermediate particles are A cl and E c baryons which 
in turn decay strongly with low Q 2 decays to A+ baryons. Furthermore, knowledge of the 
substructure gives information on the underlying mechanisms involved. 

To search for £+ + — > A+7r + and S° — > A+7r~, we require that the combination have 
a beam constrained mass within 2<jm b of the B mass peak to select events in which a B 
meson decays to a A+. We then combine this A+ with a charged pion daughter of the B 
decay and plot the M(A+7r) — M(A+) mass difference (Fig. ||). We fit these distributions 
with a Breit-Wigner function of width defined by the CLEO measurements of the E c widths 
T0[1 , convoluted with a Gaussian resolution function obtained from Monte Carlo studies, 



together with a polynomial background function. We find good evidence for both £+ + and 
TP C production in A+p7r + 7r~ and A+p7r + 7r~7r + , and for S° in A+p7r 7r~. All these signals have 
a statistical significance greater than 5a. Using analogous plots of those combinations in 
the Mb distributions outside the B mass peak, we find negligible background from true S c 
baryons that are not the daughters of the B decay mode in question. For the TP c p mode there 
are two events in the signal region which would suggest a branching fraction of the order 
of 0.25 x 10~ 4 . Our expected background in this mode, which is 0.12 events, has a 0.6% 
chance of fluctuating to the observed events. We feel that this significance, is not sufficient 
to claim to have found a signal in this mode and we prefer to present a 90% confidence level 
upper limit. There are no events consistent with the production of A cl (2593) — > A+7r + 7r~ 
or A c i(2625) — > A+7r + 7r~ in these decays which allows us to calculate 90% confidence level 
upper limits on their production. 



Table | shows the final results for the yields and branching fractions for all the modes. 
The efficiencies are calculated by our Monte Carlo simulation program. In this simulation, 
the A+ decays were generated only into the three decay modes reconstructed, using the 
measured branching ratios. To convert the quoted efficiencies to efficiencies which include 
the branching fractions of these modes, they need to be multiplied by the absolute branching 
fraction of A+ — > pK~ir + of 5.0 ± 1.3% [jO]. The yield from the S c decay modes has been 
subtracted from the non-resonant yields so that the resonant and non-resonant components 
can have different efficiency corrections applied. 

Table | includes systematic uncertainties. Major contributors to these are uncertain- 
ties due to fitting techniques, and uncertainties due to the efficiency calculation. We take 
the fitting technique uncertainty as the maximum difference obtained from different fitting 
methods. These techniques included using a scaled Mb distribution from AE sidebands for 
the background function, and fitting the AE distribution directly having first selected the 
B mass in the Mb distribution. The uncertainties from this source are 5-17%, dependent 
upon the mode. The uncertainty in the efficiency calculation is 5-8% due to uncertainties 
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TABLE I. The yields for each mode. The events in the substructure are a subset of those in 
the main modes. 



Mode 


Substructure 


Total Yield 


Substructure Yield 


A H — 




< 8 




A+pit~ 




31 ±7 


< 5.3 


A+j57r~7T + 


£ ++ p7r~ 
A+p 


110 ± 16 


14 ± 4 
23 ± 5 
< 2.3 


A + p7T _ 7T + 7r~ 
c * 


Y,®piT + 1T~ 
A+ X P7T- 


114 ± 18 


19 ±5 
12 ±4 
< 2.3 


A+pir~TT° 




76 ± 16 


13 ±4 



in the detection of the charged and neutral particles. In addition, there is a difference in 
efficiency due to possible substructure such as A*, £*,p and A intermediate states. These 
all give a slightly reduced efficiency and thus give an asymmetric systematic uncertainty. 
The systematic uncertainty due to the uncertainty in the A + — > pK~rr + branching fraction 
is expressed as a third uncertainty. 

Our limit on the branching fraction of the two body decay B° — > A+p is 0.9 x 10 -4 at the 
90% confidence level. This is tighter than the previous CLEO limit of 2.0 x 10~ 4 . A recent 
theoretical treatment by Cheng and Yang []E| using a bag model predicts a branching fraction 
of this order, whereas older theoretical predictions predicted larger numbers by a factor of 
at least 4 over our experimental limit. Our measurement of the three-body branching fraction 
B~ — > A + p7r _ and B° — > A + p7r~7r + are both consistent with, and much more accurate than, 
the previously measured modes. The three three-body decays B° — > X ++ p~7r~, B° — > £°p7r + 
and B~ — > S°p7r° have essentially identical phase-space, but only the £ ++ decay can proceed 
via both external and internal W decay diagrams, whereas the S° decays can only proceed 
via an internal W. We find the rate of all three decays to be of the same order. This implies 
that the external W decay diagram does not dominate over the internal W decay diagram, 
although naively we would expect the latter to be color-suppressed. 

In conclusion, we have measured branching fractions of B mesons into the decay modes 
A+p7r~, A+pvr~7r + , A+p7r~7r + 7r _ , and A+p7r~7r°. The first two of these confirm, with greater 
accuracy, the previous measurements. The latter two are the first observations of these decay 
modes. We find a limit on the two body decay B° — ► A+p, which discriminates between 
theoretical models. We make the first measurements of exclusive states that include S ++ or 
E° baryons. Our measurements indicate that external W diagram decays do not dominate 
over the competing internal W diagram decays in Cabibbo-favored baryonic B decays. 

We gratefully acknowledge the effort of the CESR staff in providing us with excellent 
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TABLE II. The efficiencies and branching fractions or 90% CL upper limits for each mode. 
Systematic uncertainties have been included in the upper limits. 



Mode 


Efficiency (%) 


B (10- 4 ) 


Previous Result (10" 4 ) § 


A+p 


14.9 


< 0.9 


< 2.1 


A+p7T~ 


15.8 


2.4 ± 0.6^:i? ± 0.6 


6±3 


Y? c p 


10.0 


< 0.8 




A < tp7T~7r + 


12.7 


16.7 ±1.9lJ;g± 4.3 


13 ±6 




8.0 


2.2 ±0.6 ±0.4 ±0.5 




S++P7T- 


7.8 


3.7 ±0.8 ±0.7 ±0.8 






3.2 


< 1.1 




A+p7r _ 7r + 7r _ 


9.3 


22.5 ± 2.5±£ 4 ± 5.8 


< 15 


Y,®pir + Tr~ 


5.4 


4.4 ±1.2 ±0.5 ±1.1 






5.3 


2.8 ±0.9 ±0.5 ±0.7 






1.9 


< 1.9 




A+pir~ir° 


6.8 


18.1 ± 2.9±^;g ± 4.7 


< 31 




3.8 


4.2 ±1.3 ±0.4 ±1.0 





luminosity and running conditions. M. Selen thanks the PFF program of the NSF and the 
Research Corporation, and A.H. Mahmood thanks the Texas Advanced Research Program. 
This work was supported by the National Science Foundation, and the U.S. Department of 
Energy. 
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FIG. 1. Beam-constrained mass distributions for a) A+j>, b) A+pir , c) A^pir ir + , d) 
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FIG. 2. M (A+7r)-M(A+) mass differences for combinations within 2cr of the B peak in 
the M(B) distribution, a) M(A+vr~) - M(A+) within A+pvr", b) M(A+vr-) - M(A+) within 
A+p7r~vr + , c) M(A+vr+) - M(A+) within A+pvr"7r + , d) M(A+vr~) - M(A+) within A+pvr^vr+Tr", 
e) M(A+vr + ) - Af(A+) within A+pvr^vr+vr^, f) M(A+vr") - M(A+) within A+pvr^vr . 
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